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ABSTRACT 

We present high spectral resolution VLT observations of the BAL quasar 
SDSS J0318-0600. This high quality data set allows us to extract accurate ionic 
column densities and determine an electron number density of n e =10 3 ' 3±0 ' 2 cm -3 
for the main outflow absorption component. The heavily reddened spectrum 
of SDSS J0318-0600 requires purely silicate dust with a reddening curve char- 
acteristic of predominately large grains, from which we estimate the bolometric 
luminosity. We carry out photoionization modeling to determine the total column 
density, ionization parameter and distance of the gas and find that the photion- 
ization models suggest abundances greater than solar. Due to the uncertainty in 
the location of the dust extinction, we arrive at two viable distances for the main 
ouflow component from the central source, 6 and 18 kpc, where we consider the 
6 kpc location as somewhat more physically plausable. Assuming the canonical 
global covering of 20% for the outflow and a distance of 6 kpc, our analysis yields 
a mass flux of 120 M Q yr _1 and a kinetic luminosity that is ~0.1% of the bolo- 
metric luminosity of the object. Should the dust be part of the outflow, then 
these values are ~4x larger. The large mass flux and kinetic luminosity make 
this outflow a significant contributor to AGN feedback processes. 
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Introduction 



Quasar outflows, seen as ultraviolet blueshifted absorpti on features in their spectra, 
may contribute to a variety of AGN feedback mechanisms (lElvisl 120061 ) . Such feedback 
processes are invoked to explain the: evolution of the super massive black hole (SMBH) 
rlSilk fc Reesl Il998l : iBlandford fc Begelmanl Il999h : evolution of the host galaxy (Pi Mat- 
teo et al. 2005), shape of the observe d quasar luminosity function (IHopkins et al.l l2005al ; 
Springel et al.l 120051; iMenci et al.l 120061 ) and the enrichment of the intra-cluster and inter- 
galactic media flMoll et al.ll2007n . 



The contribution of quasar outflows to AGN feedback depends mainly on their mass flux 
(M) and kinetic luminosity (Ek), both of which depend linearly on the the total hydrogen 
column density (N#) of the outflow and its characteristic distance (R) from the central source 
(see Equations 9 and 11 in § 6). Until recen tly, only loosely con strai ned estimate s have been 
available for N# and R in the same outflow. iMorris et al.l (119861 ) and lTripp et al.l (119961 ) were 
only able to find lower limits f or the distance to outflows in GC 1556+335 and HS 1946+7658, 
respectively. IWampler et al.l (119951 ) . analyzing the spectrum of QSO 0059—27 35, constrained 



the di stance of the outflow to a range spanning over an order of magnitude. Ide Kool et al. 



(120011 ) determined the distance to the outflow observed in FIRST J1044+3656 to within a 
factor of 100, and had a similar uncertainty in the inferred N# (see their table 3), which led 
to an uncertainty in M and Ek of approxi mately one ord e r of m agnitude. In the cases of 
FIRST J1214+2803 & FIRST J0840+3633, Ide Kool et all (l2002aO) were able t o constrain 
M and Ek to a similar accuracy. For the quasar 3C 191, lHamann et al.l (120011 ) estimated 
M and Ek to within a factor of a few. However, beside their internal uncertainties, all the 
works cited above suffered from unquantified systematic errors due to the use of inadequate 
absorption models that led to unreliable measurements of column densities in the observed 
troughs, which are crucial for determining almost every physical aspect of the outflows: the 
ionization equilibrium and abundances, number density, distance, mass flux, and kinetic 
luminosity. 

In order to determine accurate M and Ek in quasar outflows we launched a multi- 
stage research program: We first developed analysi s t echniques that allow us to determine 



reliable trough co lumn densities (Urav et alill999al lbl. l2002l . l2005l . boosl : IScott et al.l 12004b 



Gabel et al.ll2005l ). This was followed by identifying outflow targets in which we can find 
diagnostics for both N# and R. We searched through ~50,000 QSOs with z>0.25 and r<19.3 
in the SDSS Data Release 6 (lAdelman-McCarthy et al.ll2008l ). ~100 objects were found to 



1 Based on observations collected at the European Organisation for Astronomical Research in the Southern 
Hemisphere, Chile; under program 078.B-0433. 
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exhibit intrinsic outflow absorption systems containing troughs from excited or metastable 
levels of one or more of the following ions: He I, Sin, Fell, Fein or Nill. The key spectral 
features to obtain distances and total column densities for outflows are a) absorption troughs 
of excited and metastable levels from which we can determine the electron density, and b) 
the availability of troughs from many ions, including lines from multiple ions of the same ele- 
ment. We then began observing these objects with a combination of high spectral resolution, 
large spectral coverage and high signal-to- noise ratios t hat in agg regate allows f o r the use of 
the analysis techniques mentioned above. lArav et al.l (120081 ) and iKorista et al.l (120081 ) used 
this analysis method on VLT observations of QSO 2359-1241, deriving N# and hydrogen 
number density determinations with accuracies better than 25% and 20%, respectively. 

The object we analyze in this paper, SDSS J 0318— 0600 was found in our search through 
the SDSS data base (although first discovered by lHall et al.ll2002l ). We chose this object for a 
high resolution follow up observation, due to its high apparent brightness and its suitability 
for VLT o bservations. SPSS J0318— 0600 is a lum inous quasar (log Lb ^47.7 ergs s -1 , see §5; 
z=1.9668; lHall et al.ll2002l ; ISchneider et al.ll2003l ) that shows a myriad of outflow components 
ranging from a high ionization broad absorption line system to several narrow absorption 
line (NAL) features spanning a velocity range from —2800 to —7500 km s -1 , with respect 
to the quasar's rest frame. In a survey by Hall et al. (2002), it was classified as an iron 
low-ionization broad absorption line system (FeLoBAL) and listed as an unusual quasar due 
to the BAL located at ~27,000 km s _1 and the system of NALs, which appears as a BAL in 
higher ionization lines such as Civ. They noted that the Fell UV1, UV2 and UV3 multiplets 
were detected in absorption and that the UV38 multiplet appeared to be present as well. 



It was also noted that SDSS J0318— 0600 has a heavily reddened spectrum. iProchter et al. 



(120061 ) listed this object as a strong Mgll absorber with a re s t fram e absorption width >1.4 
A, in a search of Mgll absorbing objects. iHerbert-Fort et al.l ( 20061 ) showed that the object's 
absorption was strong, particularly in Sill A 1808, and labeled the absorber as a damped 
Lyman alpha system (DLA). However, in our VLT data, we detect ab sorption from excite d 
ionic energy levels that conclusively classify the system as an outflow (IHamann et al.lll997l ). 



The outline for this paper is as follows, in § 2, we present the VLT eschelle observa- 
tions and their reduction. For § 3, we discuss our methods of determining and measuring 
individual kinematic components. We then concentrate on a physical analysis of the main 
component of the outflow: in section 4, we present our photoionization analysis of the spec- 
trum and determine the hydrogen column density (N#), electron density (n e ), distance (R) 
and temperature (T) for the main absorption component. We derive the mass flux and 
kinetic luminosity in section 5, and in § 6 we discuss our results. 
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Table 1. VLT Observations of QSO J0318-0600 



Date 


Time a 


JD 6 


Airmass 


Seeing 


30 Sep 2006 


04:35:01 


54008.19092 


1.41 


1.08 


01 Oct 2006 


03:51:30 


54009.25268 


1.12 


1.07 


01 Oct 2006 


06:53:51 


54009.28740 


1.06 


0.97 


15 Oct 2006 


05:55:23 


54023.24679 


1.06 


1.00 


19 Oct 2006 


05:34:22 


54027.23220 


1.07 


0.82 


20 Oct 2006 


06:49:21 


54028.28427 


1.07 


0.71 


22 Oct 2006 


04:45:58 


54030.20063 


1.10 


0.66 


22 Oct 2006 


05:38:28 


54030.23505 


1.06 


0.55 



a Exposure time for each observation was 2800s 
b Julian Date presented in days +2400000 
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During October of 2006, we obtained 8 echell e spectra, each with an exposure time of 
2850 sec, of the object SDSS J0318-0600 fz=1.967. lschneider et alil2007l . RA 03:18:56.7, Dec 
—06:00:37.4) with the VLT. The observations were carried out by European Southern Obser- 
vatory (ESO) staff and we present the observational parameters in Table [TJ The spectra were 
taken using the UV- Visual Echelle Spectrograph (UVES) with a slit width of 1.0 arcsec, and 
a central wavelength position of 5800 and 3460 A for the red (Grating/Filter=CDl/HER5) 
and blue (Grating/Filter=CD3/SHP700) detectors respectively. A slit width of 1.0 arcsec 
gives a resolving power of ~40000. The data were read out in a 2x2 pixel-binned fashion, 
in high gain mode. 

The data were extracted and calibrated using the standard UVES data re duction pipeline , 



based on ECHELLE routines, which is part of the reduction package MIDAS (IBallester et al. 



20001 ) . This calibration pr ocess can introdu ce quasi-periodic oscillations in amplitude in the 
final spectral product (see lArav et al.ll2008l . on their analysis of QSO 2359-1241), which can 
create up to 2-3% systematic errors in continuum placement. The observations were cor- 
rected to heliocentric and subject to a cosmic ray removal procedure (a > 3.5). We re-binned 
the data to 7.0 km s _1 , approxi mately one r esolut ion element, and normalized the spectrum 
using the methods described in lArav et al.l (120081 ). 



2. Detection &z Measurements 



The spectrum of SDSS J0318— 0600 contains a plethora of absorption features available 
for potential diagnostics. In Figure HJ we show a sample of the VLT spectrum with line 
identifications, and present the full spectrum with identifications in our online figure. These 
identifications come from matching the Mgll doublet lines and plotting location identifiers 
for other ionic lines. We also provide an online table for all identified lines, their wavelengths, 
oscillator strengths and references. We show in Figure 1 a portion of the full online figure 
that contains readily visible lines from both the ground as well as excited levels of Fe II (lower 
energy levels up to 1873 cm -1 ). 



2.1. Identification of Kinematic Components 

We identify each of the individual kinematic absorption components in the spectrum 
using the Al II A1671 absorption line profile, which is a singlet transition with no nearby 
strong absorption features. In Figure [2], we show that there are approximately 11 separate 
kinematic components (labeled a — k), with component e clearly a blend of multiple insep- 
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Fig. 1. — Portion of the VLT spectrum showing the Fen UV 38 multiplet in SDSS J0318-0600 (the full 
spectrum can be found online). Line identifications are indicated by vertical blue lines for the expected 
position of component i, based on the redshift z= 1.9257. The redshift of the quasar is in black located at 
the top right. The log(g/) values of the identified transitions are shown above the spectrum, where g is 
the statistical weight and / is the oscillator strength. The size of the vertical blue line provides the value, 
which can be read from the y-axis scale at the top right. Each identifier lists the ion, wavelength, and lower 
energy level (in 1000s of cm -1 ). The spectrum error is plotted in solid black below the spectrum («0.05 in 
normalized flux). 



arable components. Using Al II as a template for the absorbers we search other ionic species 
for each of the 11 components. 

Next, we identify the broad, strong i component in as many lines as possible and select 
unblended troughs from which we are able to measure ionic column densities. Figure [1] shows 
the UV 38 Fell multiplet, which exhibits several absorption features from component i. We 
note that the absorption feature at ~4727 A is the k component of the Fe II A1608 resonance 
line. While this identification scheme is sufficient to locate the majority of the lines that 
show an i component, we also use the Alii as a template (see §3.3), to search for the other 
10 components in the Fell UV 1, 2 & 3 multiplet absorption complexes. 




-8000 -6000 -4000 -8000 -6000 -4000 
Velocity (km s"') 



-8000 -6000 -4000 -8000 -6000 -4000 
Velocity (km s"') 



Fig. Z. The main outflow troughs in SDSS J0318— 0600, ordered by wavelength. We identify the 11 kinematic components 

found in Alii from high to low velocities, which we label a through k. In the plots of Fell A2383, A2396, and A2600 there are 
several contaminate lines present due to the telluric A and B bands as well as blending from other Fc II troughs (identified in 
online figure). Each doublet also shows blending between components (e.g. Civ, Mgn). 



I 
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Table 2. Component i: Measured Column Densities 



Ion 


Energy Level 


Velocity Range 


Column Density 


Uncertainty" 


Number 


Methods 6 




(Ei ow cm- 2 ) 


(km s -1 ) 


(xlO 12 cm- 2 ) 


(xlO 12 cm" 2 ) 


of Lines c 




He I 


159856 


—4600 


—3716 


<90 


— 


1 


6 


Ci 





-4600 


-3716 


<37 


— 


1 


6 


Cn 


63 


-4600 


—3714 


> 19000 


— 


1 


4 


Civ 





-4600 


—3714 


29000 


3000 


1 


4 


Oi 





-4551 


—3896 


1900 


-(-zUUU 

-1200 


1 


5 


Mgi 





—4600 


—3716 


<2 





1 


6 


Mg ii 





—4600 


—3714 


3200 


400 


1 


4 


Alii 





—4477 


—3924 


400 


40 


1 


2 


Aim 





-4610 


-3833 


1560 


220 


1 


3 


Sin 





-4295 


-4106 


7220 


720 


1 


1 


Sin 


287 


-4295 


-4106 


7380 


740 


1 


1 


Si iv 





-4600 


-3714 


5600 


1300 


1 


4 


Cm 





-4295 


-4106 


14.6 


1.5 




1 


Mn ii 





-4295 


-4106 


17.5 


1.8 


1 


1 


Fei 





-4600 


-3714 


<3 






6 


Fe ii 


n 

VJ 


-4295 


-4106 


1 275 


128 


<i 


I 


Fen 


385 


-4295 


-4106 


294 


77 


3 


1 


Fen 


668 


-4295 


-4106 


64 d 




2 


1 


Fen 


863 


-4295 


-4106 


147 d 




2 


1 


Fen 


977 


-4295 


-4106 


31 d 




1 


1 


Fen 


1873 


-4295 


-4106 


163 


16 


4 


1 


Fen 


2430 


-4295 


-4106 


25 


5.4 


1 


1 


Fen 


7955 


-4295 


-4106 


8.1 


0.8 


3 


1 


Con 


3351 


-4295 


-4106 


9.4 


0.9 


1 


1 


Nin 





-4295 


-4106 


180 


18 


3 


1 


Mil 


8394 


-4295 


-4106 


64 


6.4 


3 


1 


Znll 





-4295 


-4106 


8.4 


0.8 


2 


1 



a Reflects the statistical errors from line fitting. This does not reflect the reliability of the 
measurement methods, which have increasing systematic error in order of increasing label 
number (see section 3.2). 

b Method used to determine column density, described by the respective label number; 1: 
Power Law fitting, 2: Velocity Dependent Covering Factor fitting, 3: Apparent optical depth 
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fitting, 4: "wing" fitting, 5: "tip" fitting, 6: Upper limit based on the noise method, see 
sections 3.2 and 3.3 for full details 

c Number of lines used in fitting. See Online Table for specific transition information. 

d Unreliable measurements due to weak lines 



The Fell multiplets UV 1, 2 & 3 are the best candidates to identify Fell lines from 
the other 10 components to measure, because they contain the Fell transitions with the 
strongest oscillator strengths. However, this task proves to be non-trivial in SDSS J0318- 
0600. Identifying measurable lines from a saturated and heavily blended region (see Figure 
[2]) is difficult, especially when considering the large number of strong resonance and excited 
state lines. Our strategy is to make what line identifications for all 11 components that we 
can using our Alii template and search for other lines where the blending is not so prevalent. 
We provide the full list of detected ions, transition wavelengths, their oscillator strengths 
and lower level energies in the online table. 

We find that for the lowest and highest velocity components (k and a) , enough measure- 
ments or limit determinations are available that we can effectively determine the electron 
density (n#), the ionization parameter (U#), and the total hydrogen column density (Nh, 
see §4). Unfortunately, blending prevents us from identifying enough measurable troughs for 
components b — h for proper analysis. We opt for the limiting case by taking the sum of 
these latter systems (see §3.2). 



2.2. Column Density Determinations 

With a list of measurable absorption troughs, we extract ionic column densities. We 
begin the process by measuring lines from the strongest component, (i), where we employ 
three primary methods for determining the column densities. We describe three additional 
methods for column density determination that are less reliable in section 3.3. We provide 
the measured values of ionic column densities as well as the method applied in Table [2] 

Our first method of measuring column densities is to model a set of observed Fe II reso- 
nance lines (e.g., Figure [3]) using an inhomogenous absorption model, where we approximate 
the gas distribution across our line of sight to the background source with a power law (see 



de Kool et al.ll2002d : lArav et al.ll2008l . especially Fig. 4, for a full description). The optical 



depth at a given velocity within the outflow is described by: 

T v (x) = T max (v)x a , (1) 

where x is the spatial dimension in the plane of the sky (simplified to one dimension with no 
loss of generality, see Arav et aklbood ). a is the power law distribution index, and T max {v) is 



the highest value of r at a given velocity. We simultaneously solve for both r max and a for 
each velocity resolution element of the resonance Fe II lines and apply this distribution (a) 
to the excited Fe II lines. We then apply the same distribution dictated by the Fe II lines to 
all ions because this is the simplest assumption that imposes the strongest constraints. 
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We convert the modeled r v to column density via (jSavage fe Sembachlll99ll ): 

N(v) = 3.8 x lO 14 ^-^ (cm- 2 km- l s), 
J* 



(2) 



where A is the wavelength of the line in A, / is the oscillator strength, and t v is the optical 
depth in velocity space in km s _1 , which is integrated over the spatial dimension x. Using 
Equation 2, we compute the column density for each velocity element and integrate over 
the range in velocity, which we justify later in this section. In Tabled we quote the power 
law estimation of the column densities (when available) as this is the most accurate and 
most plausibly physical explanation o f non-black saturated absorption at large distances 
( lArav et al.ll2008l ; Ide Kool et al.ll2002a ). In Tableland throughout this paper, we will refer 
to this as method 1. 

Method 2 is to model the Fell resonance lines with a velocity dependent line of sight 
(LOS) covering factor. The LOS covering fa ctor in this case is a geometric step-function 
distribution (see Figure 4 in lArav et al.ll2008l ). In other words, a fraction of the extended 
source is completely obscured by the outflow, while the rest of the source is uncovered. We 
assume that each radial velocity element has an independent LOS covering facto r, which we 



again determine using the Fe II resonance lines. We solve a fitting equation for t v ( lArav et al. 
2005h : 

I(v) = l-C(v)+C(v)e- T{v \ (3) 

where I is the normalized intensity within the trough and C(v) is the velocity dependent 
LOS covering factor. First, we solve Equation (3) for C{v) and r(v) using at least two 
resonance lines from Fell. We then calculate the velocity dependent column density using 
t(v) in Equation (2) and multiply the column density at each velocity by C(v) to obtain the 
average column density across the line of sight. Finally, we compute the total ionic column 
density by integrating across the velocity range of the trough. This method is comparable in 
accuracy to method 1 (lArav et al.ll2008l ). and the column density results are approximately 
the same as found by method 1. 

T he last of the three metho ds, method 3, is simply the apparent optical depth method 
(AOD, ISavage fe Sembachl Il99ll ). This method assumes that all photons traveling to the 
observer pass through the same amount of gas at a given velocity, which gives the equation 



for normalized intensity: I(v)=e 



-t(v) 



We again convert to column density via Equation 



(2) and integrate to find the total ionic column density. This method does not account for 
partial LOS covering and saturation effects in line centers, which can underestimate column 
densities by up to an order of magnitude. However, we find that for SDSS J0318-066, in most 
cases a LOS covering factor of 1 (i.e. the AOD approximation) yields similar N ion values as 
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those derived by the power law or partial covering method, a point we will discuss later in 
this section. 

For the upcoming photoionization modeling, we adopt the ionic column densities as de- 
termined by method 1 (power law method) with some exceptions to this due to line blending 
that are noted in Table [2j We now present our analysis for the individual components. 

2.2.1. Measurements of Component i 

We begin our column density determination by fitting the Fe II resonance lines from the 
strongest component, (i). In Figure [3], we show the fits for three Fell resonance lines using 
the power law method. We select these particular line transitions because they provide us 
with a range of oscillator strengths and are free of blending or telluric contamination. Using 
these lines, we derive velocity dependent column densities that we integrate for the total 
ionic column density. 

Unfortunately, all strong Fell lines from the lower energy level 385 cm -1 are heavily 
blended in the UV 1, 2 or 3 multiplets with lines from both the resonance and other excited 
lines (e.g., Fell* A2396, shown in Figure [2]), which renders column density measurements 
based on these transitions to have large systematic errors. However, we find relatively strong 
lines outside of the UV 1, 2 or 3 multiplets for the 385 cm -1 energy state that are useful 
such as the Fell line A1622 from the UV 38 multiplet. We also find two weaker lines from 
this energy state at A 16 19 and A2281 that are unblended. The oscillator strength values for 
these three lines are taken from Fuhr & Wiese (2006). Of the three lines, both the A1618 and 
A1622 are reported to have accurate oscillator strengths (rated as "A"g by Fuhr & Wiese), 
while the weaker A2281 has a reasonably reliable oscillator strength (rated "B"). However, we 
find that using the oscillator strengths reported by Fuhr & Wiese to be problematic. While 
the two weaker lines are well fit with these oscillator strengths, the trough depths relative 
to A1622 are unphysical. For example, scaling the A1619 trough in apparent optical depth 
by the ratio of oscillator strengths (~1.8) should be equal to the apparent optical depth in 
A 1622 (or larger in the case of partial covering), however the scaled optical depth of A 1619 
is significantly smaller than that of the A1622 trough. This means that at least one of the 
oscillator strengths could be less reliable than previously believed. Because of this issue, we 
quote an ionic column density for the 385 cm -1 energy level that is the mean of the fit to the 
A 1622 line alone and the fit to the two weaker lines with an uncertainty equal to the range 
in column density. 



2 Accuracies of "A", "B" and "C" are applied to oscillator strengths with uncertainties less than 3%, 10% 
and 25%, respectively, according to the NIST Atomic Spectra Database 
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-4400 -4300 -4200 -4100 -4000 
Velocity (km s" 1 ) 

Fig. 3. — The line fits for the three uncontaminated Fen resonance lines based on the power law method. 
The data are presented in 7 km s _1 bins with a solid line (statistical observational error bars for these data 
are also shown), while the model fit is shown by a dotted line. Each line is color coordinated with the 
wavelength located at the lower left of the plot. 



The other excited state lines from Fell (e.g., arising from 668, 863, and 977 cm -1 ) are 
also subject to heavy blending in the multiplets, but unlike the resonance and 385 excited 
state lines there is not sufficient column density to provide deep troughs with high signal-to- 
noise for the weaker lines. Thus, the column density measurements from these energy levels 
are significantly more uncertain. The exceptions to this are lines from the lower energy level 
state 1873 cm^ 1 ( 4 F) and 7955 cm" 1 ( 4 D). 

We find 5 separate absorption troughs of Fell arising from the level at 1873 cm -1 , 
which are shown in F igure HI We successfully fit these using the oscillator strengths from 



Kurucz fc Belli (119951 ) to determine the column density for this level. However, the compi- 



lation of recommended oscillator strengths of Fell from Fuhr & Wiese (2006) yields values 
lower by ~ 0.3 dex for two of the transitions, A1637.4 and A1702.0. For these transitions 
the recommended values are rated "C" , while no oscillator strengths are given for the other 
two observed transitions. In addition to the difference in the absolute values, there is also a 
difference in the relative values, and in fact our fits to the absorption troughs in the spectrum 
favors the relative values of Kurucz & Bell (1995). Thus, we opt to use the relative oscillator 
strengths from Kurucz & Bell (1995) to fit all four troughs, but the level column density 
derived this way is increased by 0.3 dex, as required by the absolute oscillator strengths of 
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-4400 -4300 -4200 -4100 -4000 

Velocity (km s" 1 ) 

Fig. 4. — The power law model fits for four lines of metastable Fen (E; ou) =1873 cm -1 ), we use the 
exponential a parameter calculated from resonance Fe II. Similar to Figure 3. 

Fuhr & Wiese. 

We find that the results from all three models (Power law, covering factor, and AOD) 
are similar for the Fell column density measurements of component i, which suggests that 
the LOS covering factor of the system is close to 1. Further evidence for a LOS covering 
factor^l case is given in Figure 0, where we show that we can use the Alii template to fit the 
Aim doublet. We scale two Alii templates (one for Aim A1855 and one for Aim A1863) in 
r space for each individual kinematic component for the blue and the red line of the doublet 
by the expected ratio of 2:1 and coadd the templates to generate a fit to the data. There 
exist only small deviations, and the overall fit to the doublet is good. The most obvious and 
important deviation is in the fit to the A1863 i component, which is due to some saturation 
in the core of the Alii template. This demonstrates that AOD is an adequate approximation 
and the LOS covering factor for all components is approximately 1.0. This same procedure 
gives similar results when applied to the Mgll doublet as well. 

2.2.2. The Physical Connection Between Material Seen in Different Troughs 

Ultimately, we want to establish the physical characteristics of the outflow based on 
the column density determinations described above. In order to do so we use the integrated 
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Fig. 5. — The Aim doublet plotted in velocity space with respect to the blue member of the doublet. 
The blue and red tickmarks and letters indicate the component positions of the Aim doublet members, 
Aim AA1855,1863 respectively. The dotted blue and red lines are the individual Alii template fits for each 
member of the doublet (corresponding to the blue and red tickmarks) . Each individual component is scaled 
in r and the resulting templates match the expected doublet ratio of 2:1 for each of the components. We 
coadd, in r space, the two Alii templates and present the resulting fit in green. The fit agrees with the data 
with the exception of the A1863 component i, which is likely due to saturation in the A1855 counterpart. 



column densities as input to photioinization models, which in turn determine the total hy- 
drogen column density and number density in the outflow. This procedure depends on two 
assumptions: 1) that the physical conditions do not change as a function of velocity across 
a given trough. This assumption allow for the use of integrated column densities instead of 
working with each velocity resolution element separately (see Arav et al 2007); 2) That the 
troughs of different ions correspond to the same outflowing gas, which is a necessary condi- 
tion for using these measurements in the same photoionization models. Both assumptions 
can be checked using the data at hand. 

The first assumption can be tested by checking the ratio measured column densities 
amongst different lower levels of the same ion vary as a function of velocity. In figure [6] we 
plot the ratios of column densities between those arising from excited levels to those from the 
ground state level of Fe II and Si II. Within the uncertainties, the ratios across the troughs 
are velocity independent, showing that the electron number density is constant across the 
flow (see § 4 for elaboration). This result gives validity to our first assumption and justify 
the use of integrated column densities for component i. We also tested this for Nill/Nill* 
and combinations between Nill, Sill and Fell resonance and excited lines and find that the 
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Fig. 6. — Component i column density ratios of Fen*/Fen and Sin*/Sin, top and bottom respectively, 
in velocity space. Within the error bars the ratios are constant across the integration range. This indicates 
that the outflowing gas has a constant electron number density across the trough, which strongly suggests 
similar ionization conditions across the velocity span of this outflow component. 



ratio is constant in all cases. This constant ratio is also seen in the results for the outflow 



system in QSO 2359-1241 as described bv lKorista et all (120081 ). 



In order to validate the second assumption we need to check the column density ratio 
across troughs from different ions and especially compare low and high ionization troughs. 
To do so we use the r(v) profile of Si II A1527 and attempt to fit the profiles of other troughs 
while using a single scale factor for each line. The results are shown in figure [71 It is evident 
that the scaled t(v) profile of Sill A1527 fits the troughs of Siiv, Alii, Aim troughs and 
Mgll remarkably well, lending strong support for the assumption that we are measuring 
the same outflowing gas. For C IV, which is the highest ionization line observed, the fit for 
the center of the trough is less conclusive since the trough is highly saturated at its core. 
However, we point out the the fit to the low velocity wing (-3950 to -3850 km s _1 ) is very 
good, which would be a coincident unless the Civ trough can be fitted with the same r(v) 
profile of Sill A1527. In summary, the data support the two main assumption that we need 
in order to use all integrated column densities in our photoionization models. 
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Fig. 7. — The top left panel shows the Sin A1527 i component. The remaining panels show the agreement 
of the Sin A1527 i component trough with the red doublet member of 5 other ions. The Sin A1527 data 
has been converted to t(v), where r(i;)=-ln(I(u)), and scaled by a single factor for each of the other lines 
in order to produce the fits (red histograms) for these troughs (data shown in black). Some discrepancy 
between the Sin A1527 line and Aim A1863 occurs near 4000 km s _1 , which is due to contamination from 
the blue k component. These agreements show that the ratio of column densities between each ion and Sin 
is constant across the profiles of the different troughs. 
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Fig. 8. — Fits to the Aim, Mgn, Civ and Siiv doublets with the composite Sin template for component 
i and the Alii template for all other components. We add two templates in a similar fashion as for Figure [5] 
These fits, however, are scaled by one constant for the blue template and scaled by 2:1 for the red template. 
We also include the Sin A1527 resonance line and A1533 excited line in the fit for C iv, which affect velocities 
greater than ~ —5500. These 4 fits clearly demonstrate that for each ion, all 11 components are present, 
and that no unexplained components appear. We note that in the C iv fit there is evidence for a component 
at —3600 km s _1 with a large optical depth whereas the optical depth of this component Alii appears to 
be significantly lower, thus the template fit is poor. The fits would improve by scaling each component by a 
different factor in optical depth. 

To demonstrate the tie between the high and low ionization troughs for the other outflow 
components, we show in Figure [8] that the combination of the Sill composite template 
(described in section 3.3 for component i) and the Alii template (for all other components) 
can be used to recreate the doublet troughs of Mgll, Aim, Siiv and Civ. Similar to 
Figure O we add two templates together in r(v). We scale the entire blue template (all 11 
components) by one factor and scale the respective red template by 2:1 in t(v). This method 
gives the fewest possible constraints for the fit. We find that all 11 kinematic components 
are present in both low and high ionization troughs, and that there are no unexplained 
kinematic components in the higher ionization state lines, which is indication that this is the 



same outflowing gas for both the high and low ionization troughs. We also note that these 
fits could be better if we scaled each individual component for the blue member and kept 
the 2:1 ratio to the red, as we did for the Aim fit in Figure El 

2.2.3. Measurements of the Other Components 

The other two components with a sufficient number of measurable absorption lines are 
the k and a components. However, we are unable to utilize methods 1 or 2 due to a lack 
of unblended Fell resonance lines. For component k, only one Fell resonance line (Fell 
A1608) is measurable and unblended. Thus, we measure the column densities based on a 
LOS covering factor of 1. We further confirm this assumption of a LOS covering factor of 
1, as with the i component, by checking that the A1608 line scales well to the unblended 
portion of Fe II A2383 according to oscillator strength. Component a has no measurable Fe II 
resonance lines. However, according to the Aim doublet (Figure [5]), we find that a LOS 
covering of ~1 holds for all components and we are able to resort to the AOD method for 
our measurements. The results for components k and a are shown in Tables [3] and Table H] 
respectively. 

Due to the lack of measurable, unblended troughs for components b through h, we 
use the apparent optical depth method, and integrate the total column density across the 
combined blended troughs for components b through h, —7150 to —4600 km s _1 in velocity, 
for the Si IV, C IV, Mg II, and Al ill doublets as well as the resonance and excited lines of C II 
(A1335, A1336) and Sill (A1527, A1533). We present these results in Table E 

2.3. Estimates and Limits 

For situations where blending between components in doublets is prevalent (i.e. Civ, 
Si IV, Cn, and Cll*) or where the trough is heavily saturated (i.e., Civ), our primary 
method of column density extraction fails. Thus, we resort to other, less accurate, methods 
to estimate the column densities, which are important when constraining photoionization 
models. We develop three additional methods in order to estimate column densities in 
blended regions, which we explain here in order of decreasing accuracy. 
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Table 3. Component k: Measured Column Densities 



Element 


Energy Level 


Velocity Range 


Column Density 


Uncertainty 


Number 


Methods" 




(E tow cnT 2 ) 


(km s^ 1 ) 


(xlO 12 cm" 2 ) 


(xlO 12 cm- 2 ) 


of Lines 6 




Cn 





-3088 -2249 


1100 


200 




4 


Cn 


63 


-3088 -2249 


1950 


300 




4 


Civ 





-2893 -2641 


1297 


130 




3 


Mg II 





-2893 -2641 


192 


20 




3 


Alii 





-3047 -2494 


35 


4 




3 


Aim 





-3047 -2494 


73 


7 




3 


Sin 





-3321 -2262 


640 


150 




6 


Sin 


287 


-2893 -2641 


352 


40 




3 


Si iv 





-2893 -2641 


140 


20 




3 


Fen 





-2893 -2641 


154 


20 




3 


Fen 


668 


-2893 -2641 


9.0 


1.0 




3 


Nin 





-3199 -2502 


<120 






6 


Nil! 


8394 


-2893 -2641 


10.0 


1.0 




3 



a Method used to determine column density, described by the respective label number; 1: 
Power Law fitting, 2: Velocity Dependent Covering Factor fitting, 3: Apparent optical depth 
fitting, 4: "wing" fitting, 5: "tip" fitting, 6: Upper limit based on the noise method 

b Number of lines used in fitting. See Online Table for specific transition information. 
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Table 4. Component a: Measured Column Densities 



Element 


Energy Level 


Velocity Range 


Column Density 


Uncertainty 


Number 


Methods" 




(E low cm- 2 ) 


(km s- 1 ) 


(xlO 12 cm" 2 ) 


(xlO 12 cm- 2 ) 


of Lines'' 




Cn 





-7523 -7271 


333 


30 




3 


Cn 


63 


-7523 -7271 


577 


60 




3 


Civ 





-7558 -7243 


734 


70 




3 


Mg II 





-7523 -7271 


28.7 


2.9 




3 


Alii 





-7523 -7271 


11.6 


1.2 




3 


Aim 





-7523 -7271 


46.0 


5.0 




3 


Sin 





-7523 -7271 


101 


10 




3 


Sin 


287 


-7523 -7271 


<50 c 






5 


Si iv 





-7523 -7271 


145 


20 




3 


Fen 





-7695 -7300 


<40 






6 



a Method used to determine column density, described by the respective label number; 1: 
Power Law fitting, 2: Velocity Dependent Covering Factor fitting, 3: Apparent optical depth 
fitting, 4: "wing" fitting, 5: "tip" fitting, 6: Upper limit based on the noise method 

b Number of lines used in fitting. See Online Table for specific transition information. 

c Column limit based on the limiting ratio "tip" fit for Sill*/Sill 
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Table 5. Components b — h: Column Densities Limits 



Element 


Energy Level 


Velocity Range 


Column Density 


Number 




(E low cm" 2 ) 


(km s -1 ) 


(xlO 12 cm- 2 ) 


of Lines" 


Cii 





-6950 -4400 


>13000 


1 


Civ 





-7150 -4600 


>9900 


1 


Mgn 





-7150 -4600 


>910 


1 


Alii 





-7150 -4600 


>440 


1 


Aim 





-7150 -4600 


>840 


1 


Sin 





-7150 -4600 


>3200 


1 


Si iv 





-7150 -4600 


>1900 


1 


Fen 





-7150 -4600 


>600 


1 



a Number of lines used in fitting. See Online Table for specific transition information. 
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In the first of our estimation methods (4), we opt to employ a template to map out 
approximate profiles for the troughs. We utilize two templates for this process. The first 
template is a fit to the Alii A1671 line for all 11 components. However, due to moderate 
saturation in the core of the Alii i component, we create a composite template using the 
Sill resonance lines for the i component. 

To create this template, we use both the Sill A1808 and A1527 lines. The A1527 line 
is similar in shape to the Alii A1671 line, and while the core is moderately saturated, it 
provides us with a good fit to the wings of the trough. In contrast, the A1808 line has a 
significantly weaker oscillator strength and therefore the core of this trough is not saturated, 
however the wings are buried in the noise. We therefore create a composite template using 
the AOD column densities of the A1527 for the wings (-3820 to -3990 km s _1 and -4390 
to -4580 km s _1 ), and the AOD column densities of the A1808 line for the core (-4390 and 
-3990 km s^ 1 ). We then convert the resultant Sill column density profile to t(v) template 
for the A1527 line. 

To measure column densities, we scale this composite t(v) template to match the red 
wings of the Mgll, Si IV, and Civ doublets as well as Cll*, which are apparently free of 
blending. We assume that the column density ratio is constant across the troughs, which is 
supported by Figure 0, scale the composite template for the i component to match and then 
integrate the total column (assuming a Q os ~l) from the template. We refer to this method 
as the "wing" fitting method (method 4). The error estimation is provided by using two r 
scaling values, one underestimating and one overestimating the noise level in the trough's 
wing. We note that these measurements use only one line, the red doublet member, which 
is reflected in Tables [2] through [5j To further test the validity of this method, we used the 
composite template to fit the wing of the Sill A1527 line and we find that the template is 
able to reproduce the column density we find in our power law fitting of the A1808 line to 
within 10%. 

The next method we designate as the "tip" fitting method, or method 5, where we 
estimate the column densities in a similar way to method 4. For example, we estimate the 
column density of O I A1302 for the i component in the region containing blended lines from 
Sill A1304 and Sill* A1309. We find that the composite template fit for the i component 
agrees well with the deep trough at ~3809 A in velocity space (shown in online figure 1), 
which appears to be the tip of O I absorption. Thus, we use the composite Si II template to 
approximate the minimum and maximum amount of column density that could come from 
O I. Due to the uncertainty in continuum placement and possible saturation, we consider 
this method less accurate than method 4, which is reflected in the larger error bars. 

Our last column density estimation method (the noise method or method 6) defines 
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upper limits on the column density of ions based on the statistical error of the data. Because 



the column density of He I is important for photoionization models (see lArav et al.l 12001 



Korista et al.ll2008l ). we estimate an upper limit for the column density of the excited 2 3 s 
level of He I (He I* hereafter) in component i for He I* A3189. To do this, we scale the Sill 
template for the i component and compare it to the region where the line would be detected. 
We vary the scaling to the noise level, find the largest possible optical depth that could be 
buried without detection and integrate across the template to determine the column density. 
We also perform this analysis for Mgl A2853, Fel A2484, and Cl A1560; and quote the 
measured limits in Table [21 

Component k presents its own problems. Unlike the i component, here we effectively 
measure column densities for Fell resonance and excited levels (E=0 cm -1 to E=668 cm -1 ), 
but we lack a Sill resonance line to provide a ratio to the Sill A1309 excited line. Therefore, 
we apply method 5 with the Al II template to find limits on the Sill. Moreover, we use 
the "wing" fitting method (method 4) with the same Al II template to estimate the column 
density for C II and C II* in component k. 

In component a, the number of measurable lines is even fewer than for component k. 
The only clean, unblend ed density diagnostic is C H to C H*. As the critical density for this 



ratio is low (~10 2 cm 3 . lBfum &: Pradhanlll992l ; iFroese Fischerlll983l ). it is important to find 



another density check. Thus, using the Alii template, we estimate an upper limit to the 
column density ratio of Sill* to Sill to be 0.49 for component a. 



2.4. Systematic Uncertainty 

We recognize that there exist systematic errors that could contribute further to the 
uncertainty in our measurements. For example, our placement of the continuum, defined in 
Section 2, could result in increased errors of the Sill A1808 and Sill* A1817 column densities 
by 5%. The amount of saturation in the core of Si IV is another systematic error that could 
be a potential source of error in column density, which at maximum, could lead to an error in 
ionization parameter of 0.1 dex. Both of these are negligible compared to the uncertainties 
due to the chemical abundance choice. 
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3. Spectral diagnostics and modeling 

The observed spectrum of SDSS J0318— 0600 exhibits distinct absorption from Cn and 
C IV, O I, Mg II, Al II and Al ill, Si II and Si IV, Cr II, Fe II, Co II, Ni II, and Zn II for component 
i. In addition, we measure upper limits to the column densities of He I*, Fel, Cl and Mgl. 
As already shown, a similar absorption trough structure is shared by the strongest transitions 
of the singly ionized species, and those of intermediate ionization (Al III and Si IV), as well 
as C IV. Additionally, component i likely has the largest column density of any component 
in SDSS J0318— 0600, given that it has the largest number of absorption line identifications 
and is the deepest absorption feature for transitions of low oscillator strengths and/or low 
relative elemental abundance. Therefore, we analyze component i here and analysis of all 
other absorption components will be presented in a subsequent paper (Bautista et al. 2009, 
in prep). 

One of the main objectives of modeling the observed ionic column densities is the de- 
termination of the total hydrogen column density of the outflow Nh and its distance from 
the central source R, which are key in the determination of kinetic luminosity of the outflow 
(see Section 6). The former is obtained directly through photoionization modeling, but the 
latter is not. Instead, modeling yields the so called ionization parameter defined as: 

Uh = t 9^ = ^jl (4 ) 

ATcR z n H c n H c 

where n H is the total hydrogen density, c is the speed of light, is the surface flux of 
ionizing photons, Qh is the rate of hydrogen ionizing photons emitted by the central object, 
which is defined by: 

f°° L 

Qh = \ ~r~dv (ionizing photons s -1 ), (5) 

JlRy 

where the integration range covers the range of ionizing photons, L u is the luminosity at each 
frequency and h is the Planck constant. From this definition, Q H can be estimated directly 
from the observed luminosity (AL A of the object and some knowledge of its spectral energy 
distribution (SED), see Sections 4.2 & 4.3). Thus, the distance R can be obtained from the 
diagnosed value of Uh whenever one has an independent determination of n H , which in a 
cloud with fully ionized hydrogen is roughly equal to the electron density, n e . 

Thus, the first step in the analysis is to determine n e from its effect via collisions on the 
relative populations of excited levels observed in the spectrum. Later, we discuss the nature 
of the SED, which is to be used in the determination of Q and as input for photoionization 
modeling. After this, we will proceed to detailed photoionization modeling. 
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3.1. Electron density determination 

Under excitation equilibrium conditions the column densities of lowly excited metastable 
levels can be used as diagnostics of electron density, assuming the atomic parameters and 
excitation processes responsible for excitation are sufficiently well understood. We have 
measured column densities for the cm -1 ground state and low-lying metastable excited 
levels within Si II, Fe II, and Ni II. 

The simplest and probably most reliable of the electron density diagnostics is that from 
the ratio of the column densities of the Si II levels, as shown in Figure [H Here, the ratio 
of the excited 3s 2 3p ^3/2 a t 287 cm -1 to the 3s 2 3p 2 P 1 / 2 ground level is calculated from 



1/2 



three different sets of effective collision strengths, iDufton &: Kingston! (119911 ) . iTayall (120081 ). 
and Bautista et al. (2009, in prep.). We measure this column density ratio from the pair of 
absorption lines at 1808.0 A (resonant) and 1816.9 A (excited state). As described above, we 
have a good deal of confidence in their measured integrated optical depths, especially their 
ratio, and we will thus place the greatest weight on this electron density diagnostic. We 
caution, however, that as the weakest amongst such pairs of transitions commonly measured 
in Sill, there are moderate uncertainties in their oscillator strengths, even in their relative 
values (about 0.1 dex). Two other such pairs of transitions that lie within our spectrum, 
centered on 1307 A and 1531 A, have smaller uncertainties in their oscillator strengths but are 
saturated in our spectrum and are not useful for direct measurements. These uncertainties 
dominate over the statistical measurement errors in the column density fits to the observed 
troughs at 1808.0 A and 1816.9 A. Fi gure [9] illustrates the electron density dependence of this 
level specific column density ratio, wherein the solid and dotted horizontal lines represent 
our measurement of this ratio and an estimated range of its uncertainty, the latter dominated 
by the uncertainties in the transition oscillator strengths. We derive an electron density of 
log(n e ) = 3.3 ±0.2 (cm" 3 ). 

The various observed transitions arising from low-lying metastable levels of the Fell 
ion allow for independent electron density diagnostics. However, constraining the electron 
density from the relative populations within these metastable levels proved more difficult. 
We present these in Fig. [10] for column densities of the 385 cm -1 , 1873 cm -1 , and 7955 cm -1 
levels. These are the only levels for which we have sufficient confidence in their measured 
column densities (see Table 2) to justify th eir utility as electron density diagnostics. In blue 



are the predictions from cloudy (v06.02. iFerland et al.lll998l ) for photoionization models 



presented in Korista et al. (2008) for the main outflow system in QSO J2359— 1241 (solar 
abundances, fixed ionization parameter and total column density for which no hydrogen 
ionization front formed, and n e ~ n H within the Fe + zone). Relevant mainly for the ratios 
involving the 7955 cm -1 level, the electron density weighted average temperature within the 
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Fig. 9. — Electron density diagnostic ratio of the populations of the first excited level of Si n and the ground 
state. The theoretical line ratios are calculated at T e = 10, 000K, although the ratio is rather insensitive to 
this ch oice. The va rious curves depict the ratios obtained from the collision strengths of Dufton fc Kingston! 
(|199lL solid curve), iTayal l (jiooi, long-dashed curve), and Bautista et al. (2009, in prep., short-dashed curve). 
The horizontal lines indicate the observed ratio (heavy solid) and estimated uncertainties (dotted), which 
are dominated by the uncertainties in the oscillator strengths of the transition pair at 1808.0 A and 1816.9 

A. 



Fe + zone of these models was approxima tely 11,000 K . In re d are the predictions of the 
collisional equilibrium model Fell atom of IVerner et al.l (119991 ) for a fixed electron temper- 
ature (10,000 K) at each electron density. This difference in adopted temperature has no 
significant effect on the populations of the 385 cm -1 and 1873 cm -1 levels, while the popu- 
lation of the 7955 cm -1 level is roughly linearly dependent on electron temperature for the 
conditions of interest (e.g., ~ 10% differences are found for the two temperatures quoted). 
The horizontal solid and dashed lines are the measured values and estimated uncertainties, 
respectively. The small differences in the predictions pertaining to the 385 cm -1 level and 
the large differences pertaining to the 1873 cm -1 level are primarily due to differences in the 
radiative transition rates (A-values) for forbidden transitions used in the two sets of compu- 
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Fig. 10. — Electron density diagnostics from the relative populations of excited metastable levels of Fen. 
The ratios are: TV(Fen* 385)/TV(Fen 0), TV(Fell* 1873)/TV(Feii 0), and TV(Fen* 79 55)/TV(Feii 0). The 
variou s curves depict the predicted line ratios by pure collisional excitation according to lBautista fc Pradhan 



Verner et al 



Itai 



(|1998l red) , by full Cloudy photoionization models dominated by collisional excitation 
blue), and a model that accounts for Civ Bowen fluorescence (green). The measured column density ratios 
are depicted by horizontal solid lines, while the corresponding errors are indicated by dashed lines. We note 
that the red curve for the TV (Fen* 7955)/TV(Feii 0) ratio is located underneath the green line. 



tations. A detailed inspection of the atomic data files of the lVerner et al.l (119991 ) model atom 
in cloudy revealed that for a few dipole forbidden transitions the model adopts A-values 
whose sources are neither known to us nor cited in their publication. For all other relavent 



transi tions the Verner et al. model atom in cloudy adopts the A-values from iQuinet et al. 
Jl996h R The largest effects of the A-values of unknown origin present in cloudy are with 
respect to the population of the 1873 cm -1 level, where the differences are substantial as 
indicated in the middle panel of Fig. [10] (blue and red curves). @. 

The electron density predicted by the Sill lines (log(n e ) = 3.3 ±0.2 cm -3 , see Figure [9]) 
and that from the 385 cm -1 level in Fell (log(n e ) ~ 3.9 cm -3 ) differ somewhat in value. 
However, considering the moderate uncertainties in the relevant atomic data of the Sill and 
Fell 385 cm -1 transitions available to us as diagnostics of electron density, and the larger 
uncertainties in the measurements in the Fell 385 cm -1 troughs, we do not consider these 



3 This situation is present within cloudy versions 08 and prior, while future versions are expected to 
adopt a revised Fen model (Ferland 2009, private communication) 

4 We have verified that the effect of the present discrepancy in the atomic data for Fen on the derived 
gas density in the main outflow system of SDSS J2359— 1241 reported in Korista et al. (2008) is minimal 
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differences to be significant. 



In contrast, the predicted electron density based on the relative population of the Fell 
a 4 F 9 / 2 level at 1873 cm^ 1 , being an order of magnitude smaller, does present a possible 
conundrum. Among all of the column density measurements of transitions from excited 
states of Fe II we believe this one to be the most secure. However, the atomic data (A- values 
and collision strengths) involved in the computations of level populations of the a 4 F term 
carry greater uncertainties than those pertaining to the levels of the a 6 D (ground) term. 
This is because the a 4 F and a 6 D terms belong to different electron configurations, i.e. 
3d 7 and 3d 6 4s respectively, and the energy difference between these terms is comparable 
in magnitude to the relativistic corrections to the Hamiltonian. Thus, the computation of 
accurate transition rates among levels of these terms requires fine tuning of large atomic 
configuration expansions (see Bautista 2008 for a more detailed discussion). 

We also consider the possibility of some selective excitation process that might lead to 
a significantly lower population in Fell a 4 Fg/ 2 than would be inferred by collisions alone. 
Neither photoexcitation by the continuum radiation field nor Lya fluorescence can explain 
the observed discrepancy. Another mechanism, known to occur in spectra of circumstel- 
lar envelopes of symbiotic stars and in the enviro ns of gamma-ray bursts is Bowen flu- 



orescence of Fell by Civ emission at 1548 A (e . g. Johansson 1983 ; Eriksson et al. 2008 



Vreeswijk et al.l 120071 ; lHartman &; Johanssonll2000l ; iMichalitsianos et al.lll992l ). Pumping oc- 
curs as the 1548.19 A resonance line of Civ closely coincides with the 1548.20 A transition 
from the a 4 F 9 / 2 metastable level to the y 4 H° 1//2 level of Fell. Subsequently, this high level 
cascades down to the ground stage of the ion, so the populati on of the a f F9/2 level is 
re-distributed among other low lying metastable levels (see also IVerner et al.l 120001 ) . Fully 
self-consistent modeling of this fluorescence process requires a detailed treatment of radiative 
transfer of Civ and Fell lines within the cloud, which is beyond the scope of the present 
paper. Nevertheless, we can estimate the relative level populations of the observed levels in 
our spectra under the fluorescence scenario by artificially including its effect in the Bautista 
& Pradhan (1998) collisional equilibrium Fell model. We find that solving this discrepancy 
requires a photoexcitation rate in the 1548.20 A transition that is at least ten times greater 
than what is available through photoabsorption of radiation from the adopted dereddened 
SED (Section 4.2). We show the resulting effects on the level populations as the green 
curves in Fig. [TTJJ (i.e., the red curve becomes the green one). Although, this fluorescence 
mechanism is able to correct the population of the a 4 F 9 / 2 level to agree with the observed 
column density at the same electron density given by Si II, the pumping rate required seems 
too large even when considering the inclusion of C IV emission from the broad emission line 
region. Further investigation of this effect is left for future work. After considering all of 
these issues pertaining to the various uncertainties, we adopt the Si n-determined electron 
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density of log(n e ) = 3.3 ± 0.2 (cm 



Finally, we measure column densities from the ground 3d 9 2 D 5 / 2 level and the 3d 8 4s 
4 F 9 / 2 metastable level at 8394 cm -1 for the Nill ion. We find that the population of the 
metastable level exceeds the predictions of pure collisional excitation conditions for any 
value of the electron de nsity. This is to be expected bec ause Ni II is very efficiently excited 
by continuum radiation (jBautista et al.lll996l : lLucylll995l ). As a consequence Nill cannot be 
used as density diagnostic here, but instead serves as a check for the flux of photons below 
the hydrogen ionization threshold from the central source. As such, Ni II can be used as a 
direct distance indicator to the outflow (see Section 4.4). 



3.2. The Spectral Energy Distribution 



The SED incident on the outflowing gas has important consequences for the ionization 
and thermal structures within the gas. In what follows, we provide motivation for our choice 
of two standard SEDs to use in the photoi onization models. While AGN are routinely 
observed in the optical/UV and X-rays (see lElvis et al.l 11994 ). the critical EUV spectral 
region (~1 Ry to 0.1 keV, rest frame) is largely inaccessible except near its low and high 
energy extremes. The present working model is that a thermal accretion disk spectrum 
dominates the energy contained within the cont inuum spectrum feature kn own as the "UV 
bump", sometimes called the "big blue bump" (jShieldslll978l ; lMalkanlll983l ). Such a feature 
is expected to peak (in units of vFy) at energies of a few Ryd berg, depending on the mas s 
of the supermassive black hole and the mass accretion rate (IShakura fc Syunyaevl Il973l ) . 
Mathews fc Ferlandl (119871 ) had this in mind when they constructed a simplified piece-wise 
power-law representation of a typical QSO SED, basing the relative strength and roll-over 
energy of the UV-FUV bump on measurements of the strengths of the He II broad emission 
lines (photon counting) and overall energy balance within the broad emission line gas (e.g., 
comparing the strength of Lya 1216 to the strengths of strong cooling lines such as Civ 
A1549), and joining the turn-down onto a power-law representation of the X-ray spectrum 
for typical X-ray bright quasars (samples of which at that time were dominated by radio 
loud QSOs). This spectrum, hereafter referred to as MF87, has since been the standard SED 
in modeling photoionized gas associated with AGN. We adopt the MF87 SED as one of two 
standards for our photoionization models. 

By the late 1990s, improved UV and X-ray observations of AGN led several investigators 
to reconsider the importance of accretion disk emission to t he FU V spectra of AGN. In 
a sample of luminous QSOs with HST spectra, IZheng et al. Jl997h foun d surprising l y sof t 
continuum SEDs at wavelengths shortward of ~1000 A. Shortly thereafter. lLaor et al.l ( 19971 ) 
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reported that the soft X-ray spectra of quasars appeared to be pointing (in vF v space) at the 



soft FUV s 
infe rred by 



opes found by Zheng et al . and proposed that a large fraction of the UV bump 



Mathews fc Ferlandl (119871 ) and standard accretion disk models was missing (see 



also lElvis et al.lll994l ). This led to investigations of cooler accretion disk models with thermal 
spectra that are Comptonized by a ~ 10 8 



(120021 ) followed up on the work of Zheng et al. with a larger sample of objects with HST 



Blaes et al. 


2001 


)• 


Telfer et al. 



spectra, which improved the statistics at the shorter wavelengths. They found similarly but 
not quite as soft FUV slopes, and concluded for radio quiet quasars "that it is plausible 
to represent the entire typical ionizing co ntinuum from ~ 10 eV to ~ 2 keV by a single 
power law" (a ~ —1.6). With this in mind lHamann et al.l (120011 . |2002| ) constructed a simple 
piece-wise power law radio quiet quasar SED to model intrinsic AGN absorption systems 
and broad emission lines. We adopt a slightly modified version of their SED as our second 
standard: from 1 micron to 1 Ry, a = —0.6; from 1 Ry to 2 keV, a = —1.6; from 2 keV 
to 50 keV, a = —0.9 and a high energy break beyond. This SED has an a ox ~ —1.43, 
very nearly equal to that of the MF87 SED (a ox _ ~ —1.40), while effectively excluding the 
presence of a UV bump peaking at FUV energiegj We refer to this as the " UV-soft SED" , 
and we show it along with the MF87 SED in Figure [TU 

With regards to these observational studies of quasar SEDs, the reader should bear in 
mind that significant uncertainties remain in the corrections in the UV for Galactic extinction 
(primarily in the adopted value of Ry), in the corrections for atomic absorption within the 
intervening IGM (mainly the Lya forest), as well as in the quasar host galaxy extinction. 
No attempt was made to correct for the last of these in any of the above works. Thus from 
an observational point of view, the existence and magnitude of the 1000 A brea k and the 
apparen t roll-down toward the observed X-rays, as suggested in lLaor et al.l (119971 ). remains 
unclear. IScott et al.l (j2004al ). in a sample of FUSE spectra of luminous Seyfertl galaxies and 
low redshift quasars, reported no FUV spectral break down to ~600 A. However, some of the 
differences between this finding and that in the more luminous AGN sample of Telfer et al. 
(2 002) may be intrinsi c, du e to differences i n the accretion disk physics. We refer the reader 



to iShang et al.l (120051 ) and iBonning et al.l (120071 ) for recent detailed comparisons between 
observations and modern accretion disk model predictions, as well as critical evaluations of 
the current uncertainties. We consider it plausible that actual quasar FUV/EUV SEDs fall 
somewhere between our two adopted standards. 



5 Hamann et al. 2001 terminated this power law portion at 1 keV, resulting in a slightly flatter alpha Q 
= -1.35. 



-32 - 



0.5 







00 

■2-0.5 



-1.5 



HISiII CII SiIV CIV 



UV soft SED 



MF87 SED 



1 

log 10 (hi/) (Ry) 



Fig. 11. — The incident SEDs considered in photoionization modeling. Both SEDs shown have the same 
bolometric flux. Here we also mark the ionization energies of some ions of interest (i.e. Hi, Sin, Siiv, Cn, 
Civ). 



3.3. Reddening of the SED 

In the previous section, we constructed a plausible UV-soft SED for SDSS J0318— 0600. 
However, it is known that the UV spe ctrum of SDSS J0318— 0600 is strongly affected by 
extinction due to dust. lHall et al.l (120021 ) found that if the Small Magellanic Cloud reddening 
curve, typical for QSOs, were adopted the extinction correction E(B — V) would be between 
~0.1 and ~0.4 depending on whether the 2000-3000 A or 1250-2000 A rest frame wavelength 
ranges were fit to the SDSS composite QSO spectrum. Because the location of the bulk of 
the reddening of the SED is unknown, we must consider a scenario in which the SED incident 
on the outflow has been reddened by dust grains. 



Figure [121 shows a portion of the SED centered around the rest frame energy range of our 
observations. Here we show the unattenuated UV-soft SED de scribed in the previous s ection 
and the FIRST Bright Quasar Survey composite spectrum of iBrotherton et al.l (120011 ) . The 
observed spectrum of SDSS J0318— 0600, also shown in the figure, is clearly affected by 
extinction. Five reddened versions of the UV-soft SED are also shown here; two of them are 



a result of the standard extinction curve of the Small Magellanic Cloud (jHutchingsl Il982l ) 
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with E(B — V)=0.1 and 0.4, and additional curves generated from cloudy models of pure 
silicate grains with a distribution of the form 

dn gr = Cn H ( — ) da, 20A < a < 0.25/xm. (6) 

These models were simulated with cloudy as the transmitted spectrum of a cloud with a 
Galactic dust to gas ratio, log n H (cm~ 3 ) = 4.0, log N H (cm~ 2 ) = 21.8 and log U H = +2.0. This 
high ionization parameter was chosen in order to minimize atomic features in the resulting 
spectrum and to allow grain extinction to dominate the attenuation. A cloud of this number 
density and ionization parameter would be located ~ 30 pc from the central source roughly 
with the expected realm of the putative clumpy, dusty torus. The dust temperature in this 
cloud has a maximum value of 400K, which is below the s ublimation temperatu re of silicate 



grains (~ 1400 K for silicates). See iPier fc Voitl (119951 ); iNenkova et al.l (120081 ) for a more 



extensive discussion of the properties of illuminated dusty tori associated with AGN. No 
graphite grains are included in the models because they lead to the 2175 A peak, which 
is absent in the spectrum. On the other hand, a large absorption peak from silicates at 
approximately 690A (~1.3 Ry) does correlate with the observed spectrum. The depth of 
such a peak is dependent on the size distribution of grains (i.e. in the value of a). A 
problem arises in determining the values of a and total extinction that are proportional to 
C, because from the UV spectrum alone no unique solution for both values can be found 
simultaneously. Fortunately, SDSS J0318— 0600 was also observed by the Two Micron All 
Sky Survey (2MASS). The photometric measurements in the J, H, and K bands (rest frame 
4162, 5602, and 7277 A respectively) are marked in Fig. [T21 It is important to note here 
that the flux in the J band may be affected by B.5 and other emission lines in the rest frame 
of the quasar. With this large wavelength range we are able to determine the magnitude 
of extinction, while the value of a is determined by fitting the slope of the SDSS spectrum. 
We find EiB — V) = 0.15 for an extinction curve with a ratio of selective to total extinction 
R v = A V /E(B -V)= 4.4 and a = 3.4. 

The extinction curve of SDSS J0318— 0600 with Ry = 4.4 is quite different from that of 
the SMC with Ry ~ 2 and the diffuse Galactic Interstellar Medium (ISM) with Ry = 3.1. 
The steepness of the present curve is more similar to that of dense molecular clouds and star 
forming regions with Ry « 5. These large values of Ry are characteristic of a distribution of 
grains that favors larger grains than in typical ISM. This would result if either large grains 
are allowed to accrete in a dense environment shielded from UV radiation or small grains 
are efficiently destroyed in a violent event, e.g. shocks. Such may be the environment of the 
dusty clouds composing the clumpy tori of AGN. 
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Fig. 12. — Spectrum of SDSS J0318-0600. The observed SDSS spectrum is depicted by the red curve and 
the 2MASS J, H, and K photometric measurements are indicated by red crosses. The black curve depicts 
our adopted UV-soft SE Dand the blue curve shows the FIRST Bright Quasar Survey composite spectrum 
(jBrotherton et al.ll200lh . The green lines show the reddened UV-soft SED by our distributed silicate grain 
model with a = 3.4 (solid line), 3.3 (dotted line), and 3.5 (dashed line). The cyan curves depict the UV-soft 
SED reddened by the SMC extinction curve with E(B — V) = 0.1 (dashed line) and 0.4 (dotted line) 



3.4. Bolometric luminosity and number of ionizing photons of 

SDSS J0318-0600 

Once the extinction that affects the observed spectrum has been determined we are in 
the position to compute the bolometric luminosity of the QSO. To do that we adopt the 
UV-soft SED shown in Figure [12] of the previous section and integrate over the whole energy 
space. This gives L Bo i = 5.0 x 10 47 ergs s _1 , assuming cosmological parameters of H Q =70 
km s" 1 Mpc" 1 , fi m =0.3, and fi A =0.7. 

With the above normalization and our adopted (UV-soft) SED we find a hydrogen 
ionizing photon rate of Q = 3.6 x 10 57 photons s _1 . However, it is unclear at this time 
where most of the extinction occurs with respect to the outflow. If the grain extinction of 
the SED occurs interior to the outflow, Q should be determined from our preferred dust 
reddened model (a = 3.4 and E(B — V) = 0.15). For this case, the effective Q — 2.1 x 
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10 56 photons s _1 . We will later use Q in conjunction with the derived values of gas number 
density and ionization parameter to estimate the distance of the outflow from the central 
source. 



3.5. Photoionization modeling 

We start by assuming constant total hydrogen density "clouds" of solar abundances 
from Lodders (2003). The optimization of all our models is based on the column densities 
of Si II and Si IV. This is because they have constrained column density measurements from 
the spectrum, expand the whole range of ionization stages observed, and they come from 
the same chemical element, thus their relative strengths are independent of the chemical 
composition of the cloud. 

In trying to reproduce the measured ionic column densities we have constructed a num- 
ber of photoionization models exploring the physical conditions and our assumptions. For a 
qualitative description of the ionization structure of the cloud see Korista et al. (2008). It 
is important in all models that the average electron density predicted by the model coincide 
with that diagnosed from the Sill lines in §4.1 within 0.1 dex. For every model, we compute 
the electron density averaged over the Si II trough, as: 



I N H N Sin n e {N H ) d N H 



< n e > Sl n= rNrr — . (7) 



This value is presented in Table |6] together with the adopted n H . 

A summary of the results of different models is presented in Table |6j Here the first 
three column densities indicate the species observed, their measured column densities and 
respective uncertainties. The subsequent column densities present the differences between 
observations and the results of various models as described below. Here we also present the 
optimal Uh and Nh determined for every model. 

We first try to model the cloud using solar abundances and two possible SEDs: our UV- 
soft SED described above and the UV-soft SED reddened as in the previous section. While 
the measured silicon column densities can be matched in both cases, this only happens when 
going very deep into the ionization front where n e and temperature have started to drop. 
Therefore, we see that in order to reproduce a log(< n e >) = 3.3 in the Sill region the 
particle density of the cloud must be log(< n H >) = 3.8. Another concern about these 
models is their predictions for O I column densities, which come to be about two orders of 
magnitude greater than allowed by observations. These results are shown as models A and B 
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in Table [6] for the unreddened and reddened SEDs respectively. Because I is tightly driven 
by Hi through charge exchange, the observed column density of O I is a direct indicator of 
the depth of the ionization front with respect to the total column density of the cloud. The 
same argument is true for the column density of Fe II. The problems with this model are 
illustrated in Figure [131 Despite the fact that our measurement of O I is poor the discrepancy 
between this and models with solar composition is significant. The current measured column 
density is 1.9 xlO 15 cm -2 with an uncertainty of a factor of 2—3. On the other hand, models 
with solar composition overestimate the column density by 1.9 dex, i.e. 1.5 xlO 17 cm" 2 . This 
would be the largest column density of all ions in the spectrum, about 5 times that of C IV. 
Despite the oscillator strength of O I being smaller than that of C IV, a column density as 
predicted would yield a very strong and completely saturated trough in the spectrum. This 
is clearly not the case. The fact that the present model goes too deep into the ionization 
front to reach the high observed ratio of Sill to Si IV column densities indicates that the 
abundance of Si relative to O must be higher than solar. 

Similar to Si, the Fe/O and Ni/O abundance ratios must also be supersolar because 
the ionic fractions of Fell, Nill, and Ol are tied to neutral hydrogen through charge ex- 
change. Thus, the Fell/Ol and Nill/Ol ratios are nearly independent of the physical con- 
dition of the cloud and can only be varied by changing the abundances of these elements 
(IBautista fe Pradhanlll998l ; iBautista et al.lll996l ). Furthermore, we find that for solar abun- 
dances of all other elements, the minimum enhancement of Si, Fe, and Ni in gas phase that 
satisfies the are [Si/H]=+0.8 and [Fe/H] = [Ni/H]=+0.3. 

We note that this abundances pa ttern is similar to that expected in high metalicity 
bulges of galaxies. iBallero et al.l (120081 ) studied the evolution of chemical abundances of spi- 
ral bulges hosting Seyfert nuclei by modeling the efficiency of star formation. They predict 
the abundances of various elements for Zj Z & = 4.23, 6.11, and 7.22. We present these abun- 
dances, which we use for photoionization modeling in Table We present these elemental 
abundances for the first and last of the metallicities, and use them in the photoion i zation 
modeling to investigate the effects of a realistic span in gas metallicity. IBallero et al.l (120081 ) 
do not predict abundances for Al, which is present in our spectra, thus we simply scale 
this linearly with Z to [Al/H]=+0.50 and 0.86 at Z/Z Q = 4.23 and 7.22, respectively. We 
use the [He/Hi vs. Z rel ation from those predicted in the chemical evolution models of 
(IHamann fe Ferlandl 119931 ). The remaining elements through Zn untabulated by Ballero et 
al. were assigned abundances as estimated by Z. 
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Table 6. Measured and model predicted absorption column densities for the outflow 

component i in QSO J0318-0600 





log(AT) (cm -2 ) 


Alog(iV) 


A 


B 


c 


jj 


E 


p 


Iog(£/ff) 






-2.25 


-2.47 


-2.63 


-3.02 


-2.85 


-3.13 


log(N H ) 






20.96 


20.92 


20.20 


20.07 


19.86 


19.73 


Iog(n H ) 






3.8 


3.7 


3.3 


3.4 


3.3 


3.3 


< log(n e ) > a 






3.3 


3.3 


3.3 


3.3 


3.3 


3.3 


< logT >sui 






4.0 


4.1 


3.8 


4.0 


3.9 


3.9 


Iog(JV H l) 






20.50 


20.43 


18.50 


19.36 


17.48 


18.26 


Ion 


Observed 




log(iV x ) f 


redicted 


- log(JVx 


^Observed 






Iog(iVx) 


Alog(iV x ) 














Si II 


16.16 


0.15 b 


+0.1 


0.00 


-0.01 


0.00 


0.00 


0.00 


Si IV 


15.75 


0.11 


+0.02 


0.00 


0.00 


0.00 


+0.01 


+0.01 


He I* 


<14.0 




-0.1 


-0.4 


-0.3 


-1.0 


-0.9 


-1.3 


Cn 


>16.3 




+0.6 


+0.5 


-0.3 


-0.1 


-0.5 


-0.5 


Civ 


16.46 


0.05 


-0.15 


-0.12 


-1.18 


-1.37 


-1.43 


-1.62 


Oi 


15.3 


+0.5 
-0.3 


+1.9 


+1.8 


0.0 


+0.9 


-0.9 


+0.3 


Mgi 


<12.3 




+1.5 


+1.3 


+1.1 


+1.0 


+1.0 


+1.0 


Mgn 


15.51 


0.06 


+0.25 


-0.39 


-0.21 


-0.84 


-0.29 


-0.66 


Alii 


14.60 


0.04 


+0.07 


-0.73 


-0.10 


-0.89 


-0.10 


-0.67 


Aim 


15.19 


0.06 


-0.60 


-1.05 


-0.73 


-1.15 


-0.65 


-1.00 


Fcii c 


15.38 


0.10 


+0.65 


+0.67 


-0.04 


+0.54 


-0.69 


+0.21 


Nin c 


14.36 


0.30 b 


+0.45 


+0.49 


+0.04 


+0.46 


-0.57 


+0.23 


Zn II 


12.92 


0.02 


-0.56 


-1.59 


-1.26 


-1.39 


-1.48 


-1.67 



Model A: UV-soft SED and solar abundances; 
Model B: Reddened UV-so ft SED and solar abu ndances; 
Model C: UV-soft SED and lBallero et all ||2008|) abundanc es for Z/Z® 



Model D: Reddened UV-soft SED and 



■■ 4.2, see Table [7] 



__.Ballero et al.l fepOSl) abundances for ZlZa, = 4.2. see TablelTl 
Model E: UV-soft SED and lBallero et all (bOOSl) abundanc es for Z/Z & = 7.2, see Table [7] 
Model F: Reddened UV-soft SED and lBallero et al.l ll2008h abundances for Z/Zq = 7.2, see Table [7] 



a Determined from the population ratio of Sill excited (287 cm -1 ) to ground levels 
b Column density error dominated by the uncertainty in oscillator strength 

c Calculated from theoretical population balance models from the column densities of observed levels 
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Fig. 13. — Column densities, temperature and electron densities in models of component i of the outflow 
of QSO J0318— 0600. In the upper panels we show in solid lines the ratios of predicted total densities of 
some species to the measured column densities vs. N H . The dashed lines are neutral hydrogen fractions. 
The results of two different models arc presented, model B, which uses solar abundances, a reddened SED 
and log Uh—-2A7, (left panels), and model F, with the Ballero et al. Z/Z Q — 7.2 abundances, also with a 
reddened SED and log J7ff=-3.13, (right panels). It is found that when solar abundances are adopted the 
measured Si column densities require a total column density that extends beyond the ionization front, which 
then results in overestimated O I and Fe n column densities 
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Table 7. Abundances used for 
Photoionization Modeling 



Ion 


Z/Z =4.2 


Z/Z Q =7.2 


[Fe/H] 


+0.83 


+ 1.07 


[O/H] 


+0.18 


+0.48 


[Mg/H] 


+0.38 


+0.69 


[Si/H] 


+0.79 


+ 1.13 


[Ca/H] 


+0.63 


+ 1.06 


[C/H] 


+0.15 


+0.30 


[N/H] 


+0.28 


+0.87 
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The results of these models are given as models C through F in Table [61 It can be 
seen that these abundances yield better agreement with respect to the measured column 
densities than when solar abundances are used. Moreover, the abundance enhancements for 
Z/Zq = 4.23 still yield overestimated column densities for O I by one to two sigma, depending 
on the SED used. On the other hand, the models with abundances for Z/Z Q = 7.2 agree well 
in terms of the O I predictions if the reddened UV-soft SED is adopted, but underestimate 
the O I column density when the non-extinguished UV-soft SED is used. 

Given our arguments for metalicities greater than solar in the outflow and the fact that 
the solar models are unable to match the n H , the true metalicity in QSO J0318— 0600 is 
likely to be several times solar. Regardless of this we are able to well constrain U H and 
N H . Models C and E constrain log(E/ H ) to -2.75 and -2.85 and log(i\y to 20.16 and 19.86 
(cm -2 ), respectively, for the unreddened SED. While, models D and F constrain \og(U H ) 
to —3.02 and —3.13 and \og(N H ) to 20.07 and 19.73 (cm -2 ), respectively, for the reddened 
SED. 

One significant discrepancy between all models and observations is in the carbon column 
densities. As indicated in Table 7, the [Si/C] ratio is rather large in these two elemental 
abundance models, and a relatively larger carbon abundance would reduce the discrepancies. 
Ballero et al. point out that their prediction for C is likely uncertain as it depends sensitively 
on the mass of the bulge. However, it is important to point out that our measurements for 
carbon are problematic, as we only have a lower limit for the C II column density and large 
portions of the C IV troughs are highly saturated. 

A relatively small uncertainty in Uh and Nh arises from the adopted standard SED. 
For example, should we adopt the MF87 SED rather than the UV-soft SED constructed in 
Section 4.2, our results for Uh and N H would be affected at a ~0.1 level. However, this 
effect is very small compared to the uncertainties resulting from the dust extinction of the 
SED and the adopted chemical abundances. 

Another discrepancy in the models stems from the predictions for Al. Most of our models 
underestimate the column densities of Alii and Aim. The predicted column density ratios 
of Al II to Al ill differ from observations by approximately a factor of four. Again, it could 
be that the true Al abundance is higher than assumed, which is not surprising as there are 
only very rough approximations on how this element scales with metalicity. But, beyond the 
absolute column densities there may also be differences in the predicted ionization structure 
of Al. The other problem is that the predicted column densities of Mg I always exceed the 
upper levels derived from our spectra. It is clear that the ionization structure of neutrals 
that form by recombination across the ionization front is not well described by our models. 
The cause for this problem is a combination of inaccuracies in the atomic parameters for 
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these species and limitations of photoionization modeling codes to treat neutral ions with 
low ionization potential (<10.2 eV). The exception to this is He I*, which is consistent with 
the models predictions. This problem is intrinsic to all modeling codes that treat radiative 
transfer through the escape probability approximation, as neutral species with low ionization 
potential are greatly affected by Lyman a radiation. Thus, a very accurate treatment of 
radiative transfer is needed to confidently model these neutral species, which is beyond the 
scope of this paper. 



3.6. The Nill lines as independent distance indicators 

Lucy (1995) demonstrated that Nill could be efficiently photo-excited by continuum 
radiation under nebular conditions. Bautista et al. (1996) constructed a detailed 142-level 
model of the Nill system for spectral analysis, which was updated with improved atomic 
data by Bautista (2004). 

Figure [T41 shows the level populations ratio between the excited level at 8394 cm -1 and 
the ground level. The upper panel shows the ratio vs. n e under pure collisional excitation 
conditions. It is noted that the uncertainty in the measured ratio column densities is domi- 
nated by that in the oscillator strengths, particularly for troughs from the 8394 cm -1 . For 
these transitions we adopted oscillator strengths from Kurucz and Bell (1995), which in the 
case of the resonant transitions were discrepant by up to a factor of two from experimental 
determinations of Fedchak and Lawler (1999). Thus, we adopt a nominal 30% uncertainty 
in the observed ratios. Further theoretical and experimental work on the oscillator strengths 
of Ni II is warranted. 

From Figure [TH one sees that the observed ratio is consistent with the predictions of 
pure collisional excitation for a wide range in n e , that includes the conditions of the absorbed 
as diagnosed from Sill. Hence, it should be possible to set an upper limit to the continuum 
flux that also contributes to the excitation of the 8394 cm -1 level. In the lower panel of 
the figure we plot the ratio against the flux of ionizing photons $ for the same SED used in 
our photoionization models. It is important to notice that only photons less energetic than 
1 Ry participate in photoexcitation of the ion and these are unaffected by the main atomic 
opacities of the cloud. It is found that the N(Nill 8394)/N(Nill 0) ratio is strongly enhanced 
by continuum fluorescence and this drives the ratio well above the Boltzman equilibrium 
populations. Under these conditions, induced photoexcitation of highly excited odd parity 
levels is balanced by rapid spontaneous radiative decay. For much more intense radiation 
fields, induced photoexcitation exceeds spontaneous decay and the level population ratio 
reaches a plateau. It can be seen that the observed column density ratio sets upper limits of 
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log w ($/c) 



Fig. 14. — Ratio of level populations between the Nill excited level at 3d 8 4s 4 F 9 / 2 level at 8394 cm -1 
and the ground level. The ratio is plotted against log(rt e ) for pure collisional excitation (upper panel) and 
log($/c) for continuum fluorescence excitation. The observed ratio is depicted by the solid horizontal line 
and the uncertainties are shown by dotted lines and the vertical dashed lines represent the upper limits for 
log 10 (<I>/c) as constrained by the measured ratio. 

log 10 ($/c) ~ +1.5 or +0.4 depending on whether one uses a reddened or unreddened SED. 
Thus, from a total particle density of \og(n H ) = 3.3 one derives a value for the ionization 
parameter log 10 (£/#) < — 1.8 or < —2.9, for reddened and unreddened SEDs respectively. 
Both of these values are consistent with that derived through photoionization modeling. It is 
unfortunate that more stringent constrains could not be derived due to the large uncertainties 
in the oscillator strengths for the absorption lines in the spectrum. 



4. Kinetic Luminosity and Mass Flux of the Outflow 

4.1. Calculating M and Ek for the Main Outflow Component 

In the previous section, we determined the ionization parameter, number density and 
total column column density of the main outflow component, (i). Here we calculate the mass 
and energetic charateristics of the outflow. First, we determine the mass contained within 
the wind. For simplicity, we assume that the outflow is in a form of a partial thin spherical 
shell at a distance R from the central source moving radially at constant speed v. The total 
mass contained in such an outflow is: 

M = 4nfim p QR 2 N H , (8) 

where Nh is the total column density of hydrogen, m p is the mass of a proton, /x=1.4 is the 
mean molecular weight of the plasma per proton, and Q is the fraction of the shell occupied 
by the outflow. Therefore, assuming diV#/dt=0, the mass flux follows as the temporal 
derivative: 

M = 87Tfim p QRN H v, (9) 

the momentum flux is given by: 

p = Mv, (10) 

and the kinetic luminosity is given by: 

Mv 2 , 
E k = = 4iTfim p QRN H v 



(11) 



For component i: v = —4200 km s _1 and table 6 gives the Njj for each of our photoion- 
ization models. To calculate the distance R we use equation (4), where n e is determined in 
§ 4.1, Uh is taken from table 6 and Qh is given in § 4.4, log Qh (photons s _1 ) as 56.3 and 
57.6 photons s -1 for the attenuated and unattenuated SEDs, respectively. In Table M we list 
the derived values for R, and M for each photoionization model. 
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Table 8. Calculated Values for Each Photoionization Model 



Unattenuated SED Attenuated SED 



Abundances 


Solar 


4.2 a 


7.2 6 


Solar 


4.2 a 


7.2 6 


Model 


A 


C 


E 


B 


D 


F 


log U H 


-2.25 


-2.63 


-2.85 


-2.47 


-3.02 


-3.13 


log N H (cm -2 ) 


20.96 


20.20 


19.86 


20.92 


20.07 


19.73 


R (kpc) 


9.4 


14.6 


18.7 


2.9 


5.5 


6.3 


E k 10 45 ergs s~ 1 ) c 


12 


3.1 


1.8 


3.3 


0.9 


0.5 


Ek/L Bo i (fto.2) 


0.02 


0.006 


0.004 


0.007 


0.002 


0.001 


M (Q .2 M yr- 1 ) 


2080 


560 


330 


590 


160 


80 



a Ballero et al. (2008) abundances for Z/Z Q =A.2 
b Ballero et al. (2008) abundances for Z/Z & =1.2 

C Q .2 = ft/0.2 
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4.2. The Global Covering Fraction 

As shown above, we determined N H and R for the photoionization models and therefore 
each model is associated with a reliable estimates of the kinetic luminosity per unit solid 
angle. In order to obtain estimates for the total kinetic luminosity we need to address the 
remaining variable in equation (1), Q, the fraction of the full sphere occupied by the outflow. 
There is no direct way to obtain the Q of a given outflow from its spectrum as we only see 
the material along the line of sight. Statistically, Civ BALs are seen in 20% of all quasars 



(IHewett fc Foftzll2003l ). This can be interpreted as: a) all quasars have BAL winds covering 
20% of the solid angle, b) 20% of all quasars have winds covering the full solid angle, c) all 
quasars have a phase of full solid angle wind lasting 20% of their duty cycle. A combination 
of these 3 options is also possible. It is important to note that with all else equal, options 
a) and c) yield the same integrated mass flux and kinetic luminosity over the life time of an 
individual quasar, and that all 3 options yield similar integrated quantities for a sample of 
quasars. 

However, we cannot measure distances to the majority of C IV outflows, but only to those 
outflows that show troughs from excited and metastable levels. For most of our objects 
(including SDSS J0318-0600), the spectra cover such transitions only from singly ionized 
species (e.g., Fe II, Sill, Cll, and Nill). Outflows with troughs from these ions are rare 



(ITrump et al.ll2006l ). Our comprehensive search of the 50,000 brightest SDSS quasars yielded 
100 such outflows. We found that only ~ 1% of the objects that show a C IV outflow show 
troughs from these singly ionized species. We therefore have two limits to explore with 
regards to the Q of the outflow. First, outflows that exhibit troughs from these 4 singly 
ionized species excited states may represent all regular C IV outflows viewed through special 
conditions that allow for the formation of singly ionized troughs. This scenario yields Q ~ 0.2 
as is the case for the frequency of Civ BALQSOs among all quasars. The second limit is 
Q ~ 0.002, which is the minimal covering fraction needed to explain the 100 occurrences in 
50,000 objects. 

Several arguments suggest that the first scenario is closer to the actual physical situ- 
ation. Since outflows with troughs from singly ionized species are rare, we may simply be 
observationally biased by this low ionization selection effect as we lack excited and metastable 
level troughs from C IV or species with similar ionization. One physically motivated picture 
is that we are looking at a normal C IV outflow through a line of sight passing through semi- 
opaque material at the edge of the putative AGN torus. Such a scenario explains why we 
see the singly ionized species (lowering of the ionization parameter); why most of the objects 
we measure have strong dust reddening (torus dust) and why these systems are so rare (a 
solid angle that has sufficient torus material to lower the ionization parameter significantly 
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but not enough to make it a type II quasar). 

Recent observations by Aoki et al. (2009, in preparation) give support to the above 
scenario. A quasar showing strong Fell absorption troughs and strong dust reddning in 
2002, showed much less dust extinction and very little Fell absorption in 2008. Thus, the 
amount of dust along our line of sight changed drastically between two epochs seperated 
by ~3.5 years in the rest frame of the source. The interpretation is that the dust source 
has vacated the line-of-sight, therefore increasing the ionization parameter of the outflow 
and reducing the column densitiy of Fell. A 3.5 years time scale favors changes in a near-by 
obscuring material (i.e., the putative torus) over material at a kpc scale and therefore support 
both of our assertions: 1) that the dusty material seen in the spectra of these outflows is 
indeed much closer to the source than the outflows themselves; 2) that these excited-state 
singly ionized absorption outflows are naturally explained by a special line-of-sight grazing 
the putative torus. 

Further support for this scenario comes from the fact that all our objects show a Civ 
BAL, which is morphologically indistinguishable from Civ BALs in the ubiquitous high 
ionization BALQSOs. We therefore conclude that it is more likely that we are seeing normal 
VI ~ 0.2 outflows through a rare line of sight. 

A direct resolution of this issue necessitates imaging of the winds. This is a difficult task 
due to the high redshifts of the objects and the presence of a very bright quasar at the center 
of the image. Nonetheless, we hope to pursue this line of research in the future. Another 
approach is to observe at wavelengths shorwards of Lya where there are a few excited state 
troughs from higher ionization species. In particular S iv/S IV*AA1062. 66, 1072.97. This ion 
is similar to C IV in ionization structure and indeed S IV troughs are much more prevalent 
than those of the singly ionized species. Measuring N H and R for these objects will therefore 
eliminate the uncertainty connected with the rarity of outflows showing troughs from singly 
ionized species, and will yield robust measurements for M and E^. We are currently pursuing 
this project, which is challenging due to the heavy Lya forest contamination. 



4.3. The Super Massive Black Hole and Host Galaxy 

Given the bolometric luminosity and assuming that the super massive black hole (SMBH) 
is accreting near the Eddington limit we find that the amount of mass being accreted is: 

Mace = ^- = 88M Q yr-\ (12) 
e r c2 
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where we assumed an accretion efficiency (e r ) of 0.1 (jSalpeteri Il964j ) and c is the speed of 
light. This large accretion rate (the outcome of a very luminous quasar) suggest that the 
mass flux values we derive (see Table 8) are not unreasonable. We expect such an outflow to 
have roughly an order of magnitude higher mass flux than the accretion rate onto the SMBH 
since the original outflow entrained material on its journey to several kpc distance. Simple 
momentum conservation requires that an outflow which started at ~20, 000-30, 000 km s _1 to 
have 5-7 times more mass flux by the time it slowed to 4200 km s _1 . Therefore, an outflow 
outflow such as this could have began with M utfiow ~ M acc . 

Next, assuming Eddington accretion, we can derive the SMBH mass and infer a mass for 
the h ost galaxy bulge. The black hole mass derived from Eddington accretion is ((Peterson 



19970 : 



M, 



BH 



4:7iGm p c 



4.0 x 10 9 M, 



(13) 



where ot is the Thompson scattering cross-section for an electron, G is the gravitational 
constant, Mbh is the mass of the SMBH, and m p is the mass of the proton. Thus, the 
spheroidal bulge component of the host galaxy (assumed to be approximately the mass of 
the host galaxy) has a ma ss Mfe u ; 9e =3.3xl0 12 M based on a SMBH to bulge ratio of 0.0012 
( jMcLure fc Dunlopl |2002| ) . This mass can be compared with the total mass ejected by the 
outflow over the ~ 10 s year duty cycle of the quasar, which is a few times 10 10 M . 



5. Discussion 

In this section we compare the results from this object to other outflows found in the 
literature; discuss the location of the dust in SDSS J0318— 0600, which is a major cause for 
uncertainty in the derived mass flux and kinetic luminosity; elaborate on the relationship 
between the dynamical time scale of the outflow and the duty cycle of the quasar; and 
conclude by describing the possible role of these outflows in AGN feedback phenomena 

In Table [9] we compare all the analyses in the literature that can be used to obtain mass 
flux and kinetic luminosity for outflows in individual quasars. We also note that with the ex- 
ception of SDSS J0318— 0600, every absorber shows evidence for inhomogeneous line of sight 
(LOS) covering. For outflow absorption, we refer to homogeneity as the spatial distribution 
of gas across our line of sight to the quasar (see section 3.2), where a homogeneous flow would 
uniformly cover the background emission source (LOS covering factor of 1). Initially, this 
may be attributed is to the large distance of the absorber, however, the outflow in 3C 191 
shows evidence of inhomogeneous covering at a greater distance than we find component i 
of SDSS J0318-0600. 
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Table 9. Properties of Measured Outflows to Date 



Object R a log N H log U ff log E k M Reference 6 

(kpc) (cm -2 ) (ergs s _1 ) (M Q yr _1 ) 



QSO 0059-2735 


0.001 - 0.05 


>21.5 C 


-0.7 




>41.1 - 42.8 


^0-2 


1 


3C 191 


28 


20.3 


-2.8 




44.0 


310 


2 


QSO 1044+3656 


0.1 - 2.1 


20.0 - 22.0 


-1.0 - - 


6.0 


44.5 - 45.4 


74 - 530 


3 


FIRST 1214+2803 


0.001 - 0.03 


21.4 - 22.2 


-2.0 - - 


0.7 


41.6 - 43.8 


0.3 - 55 


4 


FIRST 0840+3633 


0.001 


-21.3 


<-U 


I 


>41.9 


>0.3 


5 


FIRST 0840+3633 d 


0.23 












5 


QSO 2359-1241 


3 


20.6 


-2.4 




43.7 


93 


6 


SDSS J0838+2955 


3.3 


20.8 


-1.9 




45.7 


590 


7 


SDSS J0318-0600 


6 or 17 


19.9 or 20.0 


-3.1 or - 


-2.7 


44.8 or 45.4 


120 or 450 


8 



a For relative accuracies, see Section 1. 

b l-Wampler et al. (1995), 2-Hamann et al. (2001), 3-de Kool et al. (2001), 4-de Kool et 
al. (2002a), 5-de Kool et al. (2002b), 6-Korista et al. (2008), 7-Moe et al. (2009), 8-This 
Work 

c Based on Table 5 in Wampler et al. (1995) 

d Distance derived from Fe II fluorescence and no photoionization modeling was performed 
for this object 
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SDSS J0318-0600 is also the most heavily reddened in the sample, which leaves us with 



the qu estion of the location of the obscuring dust. Based on the value from ISchlegel et al. 



(119981 ). we know that the Milky Way Galactic reddening is minimal, and thus the reddening 
is intrinsic to the quasar. There exist three possible locations of the dust with respect to the 
absorber: interior, dispersed within, and exterior. The two locations for the dust we consider 
here due to their potential impact on the modeling of the outflow are that the dust is interior 
to the absorber or the dust is associated with the absorber. If the dust is interior to the 
absorber, a plausible source of dust would be the edge of the putative obscuring torus of the 
AGN unified model (lAntonucci fc Millerlll985l ). This is supported by the object discussed by 
Aoki et al. (2009, in preperation) and furthermore would provide the rare sight line needed to 
match the frequency of detections for low ionization BAL outflows that show iron absorption 
(FeLoBAL), as discussed in section 5.1. By contrast, galactic wind imaging in nearby galaxies 
shows that dust co-exists in side the winds and may be associated with the outflows (e.g., 
M82, llchikawa et al.lll994l ). lElvis et al.l (120021 ) proposed the possibility of grains forming 



within AGN outflows, as conditions allowed. We note though that if the outflowing gas 
contains the observed dust, the amount of dust depletion must be substantially smaller than 
in the ISM or SMC in order to produce the observe ionic column densities. Therefore, while 
we cannot strongly exclude the possibility of a dusty outflow, locating the dust at the edge of 
the putative torus is physically more plausible and explains several observed features in the 
FeLoBAL population: the heavy reddening of FeLoBALs, the low reddening experienced by 
C IV BAL quasar; and the higher prevalence of high ionization BALSQSOs compared with 
FeLoBALs. 



The dynamical time for component i of the outflow is tj, = R/v 



10 e 



yr (assuming 



R = 6 kpc). Therefore, if we are se eing a representative stage of the outflow, over its ~10 8 
yr duty cycle (IHopkins et al.ll2005bl ) the quasar must eject roughly hundred times more mass 
than is seen in component i. Obvious candidates for more ejected mass are the other 10 com- 
ponents of the outflow. However, Bautista et al. (2009, in prep) show that these components 
collectively do not contribute appreciably to the total mass output for SDSS J0318— 0600. 
Earlier episodes of mass ejection may not be detectable in absorption. Let us assume that 
component i is continuing at a constant velocity and no longer entraining any significant 
mass. Within 10 7 yr, its distance will increase by a factor of ten and therefore its optical 
depth will drop by a factor of a hundred, making it undetectable in absorption. To complete 
the picture we note that in addition to the FeLoBAL system, SDSS J0318— 0600 has a high 
ionization BAL at ~27000 km s _1 . This outflow component is probably much closer and 
therefore represent a very young mass ejection episode that in time may evolve to something 
resembling comp onent i. Similar situations exist in other objects (e.g. SDSS J0838+2955; 
Moe et al.l [20Q9h . 
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Given these caveats the kinetic luminosity, such as the one we have measured in 
SDSS J0318— 0600, provides ~(0.1-0.6)f2 .2% of the bolometric luminosity as feedback. We 
note that the measured mass of component i represents only the cold phase (T=10 4 K) 
of the outflow. The hot p hase of the gas (or X-ray, T=10 6 K) has been shown in Seyfert 
galaxies (e.g., NGC 3783; iNetzer et al.l 120031 ) to contribute up to 50 x more column den- 
sity than the cold phase, therefore it is possible that the gas in the hot phase of the 
outflow in SDSS J0318— 0600 could contribute significantly to the mass flux and kinetic 
luminosity. For comparison, AGN feedback models of galaxy growth and evolution re- 
quire a mech anical energy input of ro ughly 5% of the bolometric luminosity of the quasar. 
For exa mple. iDi Matteo et al.l (120051) use this ratio to recreate the observed M^-a rela- 
tion and iHopkins et al.l (l2005af ) show in their models that the same value is required to 
'uncover' a shrouded quasar after a merger. Thus, the absorption outflow we observe in 
SDSS J0318— 0600 can be a significant contributer for AGN feedback processes. 

Beyond the effects of feedback on the host galaxy, if SDSS J0318— 0600 is representative 
of the type of quasars found in clusters, then we also see that an outflow with the kinetic 
luminosit y of SDSS J0318— 0600 wo uld contribute energy to the gas that forms the ICM. The 
models of IScannapieco fc Ohl (120041 ) show that a 5% feedback ratio would suffice to h eat the 
ICM during formation of the cluster and provide the well known L x -T relationship (IKaiser 
19861 ; lArnaud &: Evrardlll999l ) Also, related to clusters are X-ray cavitie s in the ICM, whic h 
coincide with relativistic radio jets (e.g. NGC1275 in the Perseus cluster. iFabian et al.ll2006l ). 
We suggest that the sub-relativistic outflows could be related to the bubbles or perhaps are 
the driving mechanisms for bubble inflation. This follows as the in fl ation of a bubble is a 
simpler process provided a wide angle wind ({Sternberg et al.l 120071 ) . iRafferty et al.l (120061 ) 
shows that the power needed to inflate these bubbles is ~10 45 ergs s _1 , which is similar 
to the kinetic luminosity we measure in SDSS J0318— 0600 of ~l-3xl0 45 ergs s _1 . These 
similarities in fraction of the bolometric luminosity and power output suggest that quasar 
outflows play significant roles in influencing their host galaxy and neighborhoods through 
feedback. 



-51 - 



6. Summary 

1. We present 6.2 hours of VLT/UVES high resolution spectroscopic observations of 
the BAL quasar SDSS J0318-0600. 

2. Using the Al II A1671 trough, we divide the outflow into 11 kinematic components 
ranging in velocity from -7400 to -2800 km s _1 . We utilize these components to identify the 
corresponding absorption troughs in several lines from other ions. 

3. For three of these components, we measure many ionic column densities and develop 
methods to estimate ionic column densities in blended and saturated regions. The Aim 
doublet data show that all 11 components follow the expected doublet ratio of 2:1 in optical 
depth and therefore the covering of the emission source is homogeneous. Further tests for 
components i and k with Fell resonance lines demonstrate the same result. 

4. The measured column densities of Sill and Sill*, allow us to determine a number 
density of 10 33 cm -3 for component i, which is the main and innermost component. 

5. We show that the column density ratio of Nill excited to resonance can be used to 
independently determine the flux of non-ionizing photons, which when combined with the 
assumed SED yields a direct indication of the distance. These values agree with that derived 
from Sill. 

6. We analyze the intrinsically reddened spectrum of SDSS J0318— 0600, and find that 
it requires a different extinction curve from the known Small Magellanic Cloud (SMC) curve. 
We fit the spectrum with a reddened SED due to pure silicate dust and a power law grain 
size distribution with a power of a=3A, which resulted in an Ry=AA. Based on our model of 
the intrinsically reddened SED, we estimate the bolometric luminosity (log Lb i = 47.7) and 
the number of ionizing photons illuminating the absorber with and without dust extinction 
(log Q = 56.3 or 57.6). 

7. For component i, our photoionization models suggest supersolar abundances. The 
output total hydrogen column density and the so called ionization parameter for both cases 
with and without dust extinction are N H = 10 19 - 9 or 10 2ao cm -1 and log U# = —3.1 or —2.7 
respectively. 

8. Given the gas density and ionization parameter, and an estimate of the source 
luminosity, we determine the distance to the outflow for attenuated and unattenuated SEDs 
of 5.9 or 17 kpc respectively, and determine the corresponding kinetic luminosities to be (0.7 
or 2.5) xlO 45 Q , 2 ergs s^ 1 , and the mass fluxes to be (120 or 450) Q . 2 M Q yr" 1 . We find 
that component i is the dominant source of energy and mass when compared to the other 
kinematic components, which are evaluated in Bautista et al. (2009, in prep). 
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9. The ratio E k /L Bol = {0.1-0.6)Q . 2 % for component i in SDSS J0318-0600, makes this 
outflow a significant contributor to AGN feedback mechanisms. 

We acknowledge support from NSF grant AST 0507772 and from NASA LTSA grant 
NAG5-12867. MAB received partial support for this work from the NASA Award No. 
NNX09AB99G 
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